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*c/3 I Abstract 

>-. 

(~| ■ An analytic solution for Bragg grating with linear chirp in the form of 

, C^ i confluent hypergeometric functions is analyzed in the asymptotic limit of 

long grating. Simple formulas for reflection coefficient and group delay 

are derived. The simplification makes it possible to analyze irregularities 

of the curves and suggest the ways of their suppression. It is shown that 

0> I the increase in chirp at fixed other parameters decreases the oscillations in 

Q ' the group delay, but gains the oscillations in the reflection spectrum. The 

Q . conclusions are in agreement with numerical calculations. 

VO ■ PACS 42.81.Wg; 78.66.-w 

O' 

1 : 1 INTRODUCTION 

^ . Optical filters based on fiber gratings attract particular interest because of their ap- 

plications in high-speed lightware communications [ 1 1 , fiber lasers [2 1 and sensors 
lISj. The Bragg reflector is based on periodic modulation of the refractive index 

/S ' along the line ||4l|5l|. Gratings that have a nonuniform period along their length are 

c^ ■ known as chirped. The theory of linearly chirped grating holds the central place in 

the fiber optics. Chirped grating is of importance because of its applications as a 
dispersion-correcting or compensating devices [6J. The study of linearly chirped 
grating is also helpful for approximate solution of more general problem of com- 
plex Gaussian modulation Q. The group delay as a function of wavelength is 
a linear function with additional oscillations. For applications the problem is to 
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minimize the amplitude of regular oscillations and the ripple resulting from the 
errors of manufacturing {W\ . 

The purpose of this work is to present and study a solution of the equations 
for amplitudes of coupled waves in quasi-sinusoidal grating with quadratic phase 
modulation. The solution of coupled-wave equations is derived in terms of the 
confluent hypergeometric functions. Their asymptotic expansion in terms of Euler 
F-functions makes it possible to obtain relatively simple formulas for reflectivity 
and group delay. The simplification enables analysis of irregularities of the curves 
and suggestions on the ways of their suppression. 

The paper is organized as follows. The equations for amplitudes in the grating 
with quasi-sinusoidal modulated refractive index are derived in Sec. El Their an- 
alytic solutions are obtained and compared with numerical results in Sec.|3] The 
asymptotic behavior is treated in Sec.|4| Some estimations and qualitative expla- 
nations are presented in Sec.|5l Possible methods to suppress the oscillations are 
summarized in Sec.|6l 



2 Equations for slow amplitudes 

Consider a single-mode fiber with the weakly modulated refractive index n + 
5n{z). Steady-state electric field E{z) satisfies one-dimensional Helmholtz equa- 
tion 



dz' 



+ k 



^ ^ 25niz) /6n{z) 



+ 
n \ n 



£ = 0. t = i^, (1) 

C 



where z is the coordinate, k is the wavenumber in glass outside the grating, where 
5n{z) = 0, u;, c are the frequency and speed of light. The addition to mean refrac- 
tive index may be a function with phase and amplitude modulation. A family of 
analytical solutions for amplitude modulation was obtained in 0. Below we treat 
a case of phase modulation 

^^ = 2/5cos^(^), (2) 

n 

where 6{z) is the phase, constant (3 is the modulation depth. Since /3 ^ 1 we 
neglect the quadratic term in dl]). The phase is general quadratic function 

e{z) = az'^/2 + KZ + 9o, (3) 
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Figure 1 : The statement of the scattering problem. 



where k is the frequency of spatial modulation at 2 = 0, 6*0 is the constant phase 
shift. The condition of slow phase variation is 



dz 



< K. 



(4) 



Let us introduce complex amplitudes a, b of waves running in positive and 
negative directions 



E 



ae 



ikz 



+ be 



-ikz 



Keeping only resonant terms and neglecting the parametric resonance of higher 
orders at the detuning 



q = k- k/2 <^ko = k/2, 
we get the equations for coupled waves 

a' = iko(3e-^'''^''^% y = -ikolSe^"''-''^'^ 



(5) 



(6) 



where prime denotes z-derivative. 

Set ^ conserves |ap — |6p, since the signs in right-hand sides of equations 
are different. The same equations with identical signs conserve the sum of popu- 
lations |ap + |6p and describe the amplitudes of probability in two-state quantum 
system. The exact solutions in this case are of importance in quantum optics, then 
they are studied in details ifTDlfTT l. Within the limits of resonance approximation 
^ we replace k in front of exponents © by Ajq. 

Finding the derivatives of I© with respect to z we get complex conjugated 
second-order equations 



a" — ia{z — zq) a' 



kl(3^a = 0, b" + ia{z- zq) b' - kll3% = 0. (7) 



Here zq = {2k — K)/a = 2q/a is the coordinate of resonance point for the 
wave with wavenumber k. It is the turning point where the wave with given q is 



reflected. The parametric resonance for central wavenumber q = k — k/2 = 
occurs at z = 0. Let a > 0, then for the red detuning g < we have zq < 0, in 
opposite case g > of blue detuning zq > 0. 

Consider the Bragg grating written in the interval —L ^ z ^ L. The problem 
of left reflection coefficient calculation is illustrated by Fig.[T] Boundary condi- 
tions are defined by the scattering problem statement. We set amplitude b at the 
right end equal to zero 

b{L) = (8) 

and get the reflection and transmission coefficients 

a[—L) Ci\~J^) 

The chirp is weak and satisfy Q, then the equations for complex amplitudes are 
valid when 

aL < K. (10) 

Note that set Q is symmetric under transformation a — > —a, q — »• — g, a ■^ h. 
Then the right reflection coefficient can be obtained from the expression for left 
one by changing signs of parameters a and q. 

3 Solution 

Equations Q are reduced to the confluent hypergeometric form by the substitution 

t = m{z~ZQf/2: 

ta + { t\ a + irja = 0, 

where dot denotes the derivative with respect to new variable t,r] = P'^kl/2a is the 
adiabatic parameter. The equation for second amplitude b is complex conjugated. 
The general solutions at —L < z < L are linear combinations 

a{z) = Aiui{z) + A2U2{z), b{z) = Biul{z) + B2ul{z) (11) 

of the Kummer confluent hypergeometric functions llT2ll : 

Ui= F {-IT]; i; ia{z - z^f /2) , 

U2 = {z- Zo)F (i - ir/; f ; ia(z - ZoY/2) ; (12) 

a X aia + 1) x^ 
F(a;c;.) = l + -- + -l_^- + ..., 
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where Ai, A2, Bi, B2 are constants and the asterisk denotes the complex conjuga- 
tion. The solution was obtained in Ifl^l for optical waveguide. The solution for 
coupled- wave equations with identical signs has been obtained in the context of 
nonadiabatic population inversion in two-level system [14|. 

The relations between constants can be obtained from set ® near resonance 
point z = Zq where a = Ai + A2{z — zq) + 0{z — zqY, b = Bi + B2{z — zq) + 
Oiz-zo)': 

^ = iA;o/3e'^o-^"^o/2^ ^ = ^,koPe-''°+'</'. (13) 

Bi Ax 



The right boundary condition dH) yields the ratio of coefficients A\ and A 
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Ai 1 B2 F{w,i;-m{L-zof/2) 

^ A2 kl(5^Bx f3^k^,{L-Zo)F{i + w,l;-m{L-Zoy/2y ^ ^ 

The left reflection and transmission coefficients Q can be expressed in terms of 
confluent hypergeometric functions 

^-,eo+ra4/2 ^.(_^) ^ /3XpU*2i-L) ^_ pm{L) + U2{L) ^^^^ 

ikoP pui{-L) + U2{-L) ' pui{-L) +U2{-L)' 

where functions ui 2 are defined by (fT^ and p is defined by (IT^ . 

The reflection spectrum, i.e., the reflectivity R = |rp as a function of detun- 
ing q, is shown in Fig. |2l (a). The central frequency of the spectrum comes to 
resonance at 2; = 0, in the middle of grating. The central part has a flat top at 
high adiabatic parameter, as the upper curve shows, and the maximal reflectivity 
is close to 1 . The width of central part is proportional to the length L. The reflec- 
tivity is relatively high if the turning point zq lies inside the grating \zo\ < L. This 
inequality gives the bandwidth |g| < aL/2. There is no parametric resonance at 
higher detuning, when \zq\ > L, and the reflectivity is small. Fig.|3l(a) shows how 
the bandwidth grows up with the chirp parameter a at fixed modulation depth 
p. The adiabatic parameter decreases with a, then the reflectivity in the center 
decreases from curve to curve. 

The spectrum was recalculated numerically by T-matrix approach. The ini- 
tial Helmholtz equation ([T]) was solved numerically with neither approximation of 
slow envelope, nor quadratic term {5n/nY neglecting. The number of layers per 
period of spatial modulation was fixed at A^ = 32, then the step varied along the 
grating. The spectra for n = 1.5 and the same parameters are shown in Fig.|2(b) 
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Figure 2: (a) Reflection spectrum R{q) at a = 600 cm '^,L = 0.5 cm,ko = 
6 X 10^ cm-\ from the top down (3 = /3o, f3o/2, po/A, where (3o = 0.67 x lO^^. 
(b) A part of the lower curve /3 = (3q/A, crosses denote the numerical result. 
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Figure 3: (a) Reflection spectrum from the top down at L = 0.5 cm, P = 0.33 x 
10~^, /co = 6 X 10*^ cm~^ and a = ao, 2«o, 5ao„ where «o = 600 cm"^. (b) A 
part of the lower curve a = 3q;o, crosses denote the numerical result. 



by crosses. The numerical results are very close to analytical, since both dimen- 
sionless parameters controlling the validity of coupled- wave approximation are 
small: P ~ 10"'^, aL/n ~ 5 x 10^"^. At higher parameter a the deviation of 
coupled- wave equations solutions from that of Helmholtz equation increases, but 
not dramatically, as shown in Fig.|2](b). The origin of the deviation is resonance 
approximation Q. We replace k by ko in coupled-mode equations ©, while the 
Helmholtz wave equation O involves k. Comparing Fig.|21(a) and Fig.|3](b) we 
see that the latter involves higher detuning, then the deviation is greater at higher 
q = k — ko. 

The group delay found from analytical solution (fTSt is plotted in Fig.|4l(a) at 
the same parameters as the reflection spectrum in Fig. El The deviation of curves 
from the linear dependence, the group delay ripple, manifests itself as oscillations 
with variable frequency. The frequency grows up towards the blue end of spec- 
trum in agreement with results from |fT5lj5^. For the negative chirp (or when the 
incident light enters from the right) the frequency grows up towards the red edge 
of spectrum. The maximum deviation from the averaged slope decreases with de- 
creasing modulation depth /?. Meanwhile, the ripple in reflectivity increases for 
small p. A fragment of group delay characteristics is shown in Fig. HI (b) along 
with numerical calculations. Dots obtained from numerical calculation are very 
close to the curve given by analytical formula. 

It is difficult to analyze the solution in its general form. In particular the cum- 
bersome expression for group delay, the derivative of (IT5t with respect to the 
detuning, is not presented here. Let us simplify expressions using the asymptotics 
of Kummer functions in the next section. 



4 Asymptotics 

The asymptotic expressions for the reflection coefficient can be obtained from 
(IT5t in two cases. The first case is the resonance condition at the left end, namely, 
detuning q = —aL/2 for which zq = —L. In this case it follows from (IT^ that 
Mi(— L) = 1, U2{~L) = 0, and then from (1131 

R = \r\'^ ^ taiih TCI]. (17) 

The other case is when the resonance point zq being far from both ends inside 
the grating: |g| < aL/2 and a(zo ±L)^/2 ^ 1. The asymptotic expression of the 
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Figure 4: (a) Group delay (ps) as a function of detuning q (cm^^) at a = 
600cm~^,L = 0.5 cm, /co = 6x10"^ cm"\ from the top down /5 = /3o,/9o/2,/9o/4, 
where /3o = 0.67 x 10^^. (b) Comparing with numerical calculation denoted by 
dots. 




Figure 5: The reflectivity i? as a function of adiabatic parameter r] at zq 
(solid line) and at zq = (dashed). 



confluent hypergeometric functions lIT^ at |argx| < vr 

T(c) /"e'^^'V Tie 



F{a; c; x) 



T(c — a) 



X 



+ 



T(a 



-e^x" '', 



e=<+'' '"^">°' (18) 
-1, Imx < 



allows one to simplify expression (fTSt . 

The reflection coefficient can be written using (fT2ll 



2 V a 



n -2i7? 






(19) 



where ^±(g) = a(L ± 2q/af/2, <p = arg[r(ir/)r(l/2 + ir^)], Ro = l- e"^^"^ and 
we omit terms of the order of l/aL?' -C 1. The enumerator and denominator of 
the fraction in the second line of ( IT9t are close to 1, if ^± ^ 1. Then the formula 
for reflectivity becomes simple 



R^ RQ = l-e 



-iwr) 



(20) 



Both curves dTTJ) and (l20b are shown in Fig. |5l One can see their saturation, 
moreover, when the turning point zq = —L, the saturation occurs later than when 
the turning point is far from both ends. 
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The group delay obtained from (IT^ is 



d arg r n dim In r n ( Aq 2rj 

duj c dq c\ a \ R^a 

X (l + i?o)cos(0-7r/4 + V^+ -2?7lnV^+) + 

+ (1 - i?o) cos (0 - 7r/4 + ^_ - 2r7 In ^_)] } , (21) 

where we neglect terms of the order of 1/aL^ ^ 1. Expression (|2T1) involves 
three terms. The first (in the first line) gives the averaged slope. It is a linear 
function within the bandwidth. Its slope depends on parameter a. At 1 — _Ro ^ 
1 the second term (the second line) gives the ripple, chirped oscillations. The 
frequency of these oscillations is double distance from left end of the grating to 
the reflection point zq. Their frequency ^+(q') = 2L + Aq/a grows up towards the 
blue edge of the spectrum. When reflectivity Rq becomes smaller, the last term 
(the third line in EqlSTT) proportional to Tq = 1 — -Rq comes into effect. It gives 
the additional oscillations with variable frequency il)'_ (g) = 2L — Aq/a that grows 
up towards the red edge of the spectrum. It is precisely the sum of two chirped 
oscillations with significantly different frequencies that the left part of the lower 
curve in Fig.|4|(a) displays. Magnified view of the corresponding fragment is also 
shown in Fig. |4l(b). If we change the sign of chirp parameter a, then functions 
ip± switch their roles: ^+ *-^ ip_. Therefore at high reflectivity 1 — -Rq ^ 1 the 
spatial frequency of leading oscillations ^^ = 2L + Aq/a decreases towards the 
shorter wavelengths. 

The amplitude of oscillations in group delay (1211) increases when Rq tends 
to unity, while that in the spectrum decreases. Formula for the reflection inside 
the bandwidth can be obtained from (IT^ with the accuracy to the next order of 
transparency Tq = 1 — Rq 



cos{(f) — n/A + ^+ — 2ri\mp^) cos(0 — it /A + ■?/'_ — 2ri\nilj. 



(22) 



^/^{L + 2q/a) y/aj2{L - 2q/a) 

At high reflectivity Rq -^ 1 oscillations (l22b are suppressed. The first term in 
square brackets describes oscillations with frequency 2L + Aq/a, their amplitude 
gains towards the red edge of the spectrum. The second term corresponds to oscil- 
lations with frequency 2L — Aq/a with amplitude growing towards the blue edge. 
Both approximate formulas (0^ and (OTT) for oscillations are plotted in Fig. |3 and 
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Figure 6: The group delay calculated according to asymptotic formula (12111 at 
a = 600 cm-2, L = 0.5 cm, fco = 6 x 10^ cm'^ f3 = 0.67 x 10^1 Dots denote 
the exact Kummer solution. 

Fig. El respectively. As figures illustrate, the asymptotic expressions nearly coin- 
cide with exact Kummer solutions. The departure of the simple formula from the 
Kummer solution (left edge in Fig.|6land both edges in Fig.|3) occurs when we get 
the limit of applicability of the asymptotic expansion. The turning point should 
be located far from the ends of grating, i.e., \L ± 2q/a\ ^ Lejj = (2n/\a\Y^'^. 
Asymptotics are broken when the turning point occurs too close to the end. 

The dependence on parameters a, (3 in Fig. and |3] can also be explained by 
the asymptotic expressions. At fixed chirp parameter a the adiabatic parameter 
r] = {koPy/2a in (EUt decreases with decreasing the modulation depth p. Then 
reflectivity _Ro at g = is relatively small and oscillations with amplitude 1 — -Rq in 
the spectrum become noticeable. At fixed /5, on the contrary, the adiabatic param- 
eter decreases with increasing a. It is the reason of the most evident oscillation in 
the spectrum corresponding to the higher chirp parameter a. 

5 Discussion 



The reflectivity is maximal at k = ko = k/2 ■ 
of modulation in the middle of the grating, at z 



vr/A, where A is the period 
0. The spatial frequency of 
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Figure 7: Reflection spectrum calculated by asymptotic formula (|22b at a = 
600 cm~\/3 = 0.17 X 10"^, ko = 6 x 10^ cm"^ Dots denote the exact Kum- 
mer solution. 

modulation 9'{z) = az + k depends on coordinate z. Then at some distance from 
the center the wave with k = ko comes out from the resonance. The dephasing 
occurs when 9 = az'^/2 ~ vr, i.e., at distance z = Les ~ (27r/Q;)^/^. The 
effective number of strokes along length Lgff should be large Mgff = Lgff/A = 
(27r/a;)^/^A^^ ^ 1. Moreover, to provide the high reflectivity it should satisfy 
the stricter limitation of dense grating M^ffP > 1. From here we get a condition 
for adiabatic parameter 



V 



2a 



1/4:71. 



The bandwidth of the reflection spectrum is aL, as shown in Sec.|3l The fronts 
of spectrum are determined by the effective length Lgff. When point zq = 2q/a 
is placed outside the grating at distance ~ Lgff from the end, the reflection almost 
vanishes. The width of fronts is 5q = aLgg- = 1/Lgff. The fronts are steep while 
Lgfl <^ L, i.e., in the limit of long grating. 

Phase modulation 9(z) provides the parametric resonance condition for differ- 
ent wavelengths. The shorter waves meet their resonance at longer distance zq = 
2q/a, and then the group delay of blue light is more than that of red one, Fig.|4| 
The linear dependence of the average group delay (1211) upon the detuning has also 
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-L Zq = 2q/a +L 

Figure 8: Configuration of compound cavity: left "mirror" z = —L is the left edge 
of grating, right "mirror" z = +L is the right edge of grating. Middle variable 
"mirror", the turning point z = zq, is located at different positions depending 
on the wavelength. Then the ripple frequencies are determined by the variable 
lengths of sub-cavities I and II. 

simple explanation. The delay r = //f gr is defined by double distance from start- 
ing point to the resonance for given wavenumber / ^ 2zq = 4(/c — ko)/a. Here 
Vgr is the group velocity of light. If the chirp a is negative, then the sign of delay 
characteristics becomes negative. 

The ripple outside the reflection spectrum bandwidth. Fig. 12131 with period 
tt/L are the Gibbs oscillations originated by steep boundaries, i.e., reflection from 
the grating edges. The aperiodic oscillation inside the bandwidth arise from the 
triple-mirror cavity with moving middle mirror, Fig.lU The wave reflected to the 
left from turning point zq could reflect back to the right from the left end of the 
grating. Then the cavity appears between z = —L and z = zq; its effective length 
is I = Zq + L = 2q/a + L. It results in oscillations with period 7r/(L + 2q/a). 
The cavity with variable "mirror" is longer for blue spectrum and shorter for red, 
then the frequency of oscillations increases with q, as mentioned in paper ifTSll . At 
Rq ^ I these oscillations are suppressed in the reflection spectrum, but remain in 
the group delay characteristics. If the reflectivity is not close to 1, the additional 
oscillations come into effect due to the "right" cavity with variable "left mirror". 
Their period 7r/(L — 2q/a) on the contrary is longer for red spectrum. These 
oscillations are suppressed at i?o — > 1 both in reflection spectrum and group delay 
characteristics. 



6 Conclusions 

Thus, the analysis of the reflection spectrum and group delay of linearly chirped 
grating becomes simple if the turning point zq = 2q/a is far from both ends of the 
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grating compared to the effective length Lgff = (27r/«)^/^. Formulas for reflec- 
tivity demonstrate the irregular oscillations in the reflection spectrum when the 
adiabatic parameter is not large. The oscillations are aperiodic and their ampli- 
tude slowly increases from the center of spectrum. The nature of the oscillations 
is reflection in compound cavity with a mobile middle "mirror". There are two 
terms in asymptotic expression. The first has a period vr/(L + zq) (round trip in 
the left sub-cavity), the second — n/(L — zq) (round trip in the right sub-cavity). 
The oscillations in group delay characteristics have the same origin. The differ- 
ence is that the right sub-cavity takes a negligible part in forming the oscillations 
of group delay characteristics at i?o — * 1- 

The amplitude of oscillations is suppressed at high chirp parameter a even 
at fixed reflectivity. The conservation of high reflectivity with increasing a re- 
quires increasing parameter /3. In order to suppress both oscillations one must 
choose as high the modulation depth as possible, but the limitation exists in fiber 
Bragg grating manufacturing. The alternative method to diminish the unwanted 
echo might be to provide the signal dephasing by apodization, i.e., smoothing the 
grating profile |5|. 
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